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ABSTRACT 

The correlation between black hole mass Mbh and stellar velocity dispersion 
a in nearby elliptical galaxies affords a novel way to determine Mbh in active 
galaxies. We report on measurements of a from optical spectra of 7 BL Lac 
host galaxies. The derived values of a are in the range of 160 - 290 km s _1 
corresponding to Mbh of 5 x 10 7 to 1 x 10 9 M . The average ratio of Mbh to 
the host galaxy mass is 1.4 x 10~ 3 , consistent with that estimated in other active 
and inactive galaxies. The velocity dispersions and the derived values of Mbh 
of the BL Lacs are similar to those obtained for low redshift radio galaxies, in 
good agreement with the predictions of the unified models for radio-loud active 
galaxies. 



Subject headings: BL Lacertae objects - galaxies:active - galaxies: elliptical and 
lenticular - galaxies: kinematics and dynamics - galaxies: nuclei 
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1. Introduction 

The mass of the central black hole (BH) is of paramount importance in theoretical 
models of AGN. In particular, the dependence of BH mass (Mbh) on the global host galaxy 
properties provides clues to the role of BHs in galaxy formation and evolution. Dynamical 
determination of M BH in AGN is difficult, because of the bright emission from the nucleus. 
The main method that has proved to be successful for AGN is reverberation mapping of 
broad emission lines, which is extremely time consuming and gives results on Mbh that 
depend on the assumed geometry of the accretion disk. Therefore, only for a few well 
studied quasars and Seyfert galaxies Mbh is known (see e.g. Kaspi et al. 2000; Nelson 2000; 
Wandel 2002 and references therein). Reverberation mapping cannot be employed for BL 
Lac objects because they lack prominent broad emission lines, therefore other methods need 
to be applied. The discovery of a correlation between Mbh and the luminosity of the bulge 
in nearby early-type galaxies (e.g. Magorrian et al. 1998) offered a new tool for evaluating 
Mbh (see the recent reviews by Merritt & Ferrarese 2001 [hereafter MF01]; Kormendy & 
Gebhardt 2001). This correlation has been applied so far for a sample of nearby quasars 
(McLure & Dunlop 2001) and BL Lacs (Treves et al. 2001). 

Recently, a stricter correlation was found relating Mbh with the stellar velocity dis- 
persion o of the spheroidal component in nearby inactive galaxies (Gebhardt et al. 2000; 
Ferrarese & Merritt 2000). This relationship clearly demonstrates a connection between BHs 
and bulges of galaxies and has spurred a substantial effort in theoretical modelling (e.g. Silk 
& Rees 1998; Hahnelt & Kauffmann 2000; Adams, Graff & Richstone 2001). The relation- 
ship appears to predict more accurately Mbh-, but requires the measurement of a in the 
host galaxies of AGN that is difficult to obtain, in particular for objects at moderate or high 
redshift and with very luminous nuclei. On the other hand, BL Lacs have relatively fainter 
nuclei than quasars (e.g. Falomo et al. 1999), and for them this measurement (at least for 
low redshift objects) can be secured with a single spectrum observable with a medium-sized 
telescope. 

We present here medium resolution optical spectroscopy of the host galaxies of 7 BL Lac 
objects from which we derive the stellar velocity dispersion. According to the shape of their 
spectral energy distribution (SED), BL Lacs are broadly distinguished into two types (see 
Giommi & Padovani 1995) : those whose SED peaks at near-infrared/optical and the 7-ray 
MeV regions (low frequency peaked BL Lacs or LBL), and those that have SED peaking in 
the UV/X-ray and the 7-ray TeV energies (called high frequency peaked BL Lacs or HBL). 
Our selection of nearby (z < 0.2) BL Lacs includes 5 HBL and 2 LBL. For all observed targets 
high quality images have been obtained either from the ground (Falomo & Kotilainen 1999) 
or with HST + WFPC2 (Urry et al. 2000; Falomo et al. 2000). From these images, the 
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characterization of the host galaxies and the nuclear luminosity can be obtained. 

2. Observations and data analysis 

The observations were obtained in June 2001 using the 2.5m Nordic Optical Telescope 
(NOT) equipped with ALFOSC 1 . Spectra were secured using two grisms to cover the spectral 
ranges 4800 - 5800 A (setup A) and 5700 - 8000 A (setup B) at 0.54 A pixel" 1 and 1.3 
A pixel -1 dispersion, respectively. This allows us to measure the absorption lines of H/3 
(4861 A), Mg I (5175 A), Ca E-band (5269 A), Na I (5892 A) and the TiO + Cal (6178 A), 
TiO + Fel (6266 A) and other absorption line blends from the host galaxies at a spectral 
resolution R ~3000. 

The chosen grisms combined with a 1" slit yield a spectral resolution for velocity dis- 
persion measurement of ~60 - 80 km s -1 , which is adequate for the expected range of a in 
luminous ellipticals (e.g. Djorgovski & Davis 1987; Bender, Burstein & Faber 1992) such as 
the hosts of BL Lacs. In addition to the spectra of the BL Lac hosts, we acquired spectra 
of bright stars of type G8-III to Kl-III, that exhibit a low rotational velocity (Vxsin (i) < 
20 kms -1 ). These are used as templates of zero velocity dispersion. Furthermore, spectra of 
the well studied nearby elliptical galaxy NGC 5831 were also secured in order to provide a 
test of the adopted procedure to derive a. 

During the observations, seeing ranged between 1" and 1.5". The targets were centered 
into the slit or positioned 1" away from the nucleus and then the ID spectrum was extracted 
from an aperture of 3" - 5" diameter, which is in all cases within the effective radius of 
the host galaxy. In one case (Mkn 180), spectra with the object both centered into the 
slit and off-centered by 1" were taken but no significant difference was apparent in the 
shape of the spectral features. Standard data reduction was applied to the spectra using 
the tasks available in the IRAF 2 package. The procedure includes bias subtraction, flat- 
fielding, wavelength calibration and extraction of ID spectra. For each observation, we took 
two spectra and combined them in order to remove cosmic ray hits and other occasional 
spurious signals in the detector. In Table 1 we report the list of our targets together with 
the instrumental setup and the S/N of each spectrum derived from the continuum in the 
middle of the observed spectral range. 



1 See URL http://www.not.iac.es for instrument characteristics 

2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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The stellar velocity dispersion a was determined using the Fourier Quotient method 
(e.g. Sargent et al. 1977) implemented in the IRAF STSDAS package. The spectra were 
first normalized by subtracting the continuum, converted to a logarithmic scale and then 
multiplied by a cosine bell function that apodizes 10% of the pixels at each end of the 
spectrum. Finally, the Fourier Transform of the galaxy spectra was divided by the Fourier 
Transform of template stars and a was computed from a x 2 fit with a Gaussian broadening 
function (see Bertola et al. 1984; Kuijken & Merrifield 1993 for further details on this 
method). The r.m.s. scatter of the a results using different template stars was typically 
~10 kms -1 and can be considered as the minimum uncertainty of the measurement. The 
observed values of a and their estimated errors are reported in the last column of Table 1. 

For three objects we have spectra in both spectral ranges. The resulting stellar velocity 
dispersions are in all cases in good agreement, with average difference 12 km s _1 , ensuring 
sufficient homogeneity of data taken with different grisms and/or resolution. Note, however, 
that there is a tendency for the red grism (lower resolution) data to result in slightly larger 
value of a. For the nearby elliptical NGC 5831 we obtained a = 167±5 and 185±10 km 
s _1 for the setup A and B, respectively. These values are in good agreement with previous 
measurements in the literature (< a > = 168 km s -1 ; Prugniel et al. 1998). In Fig. 1 we 
show an example of the spectrum of the BL Lac object Mrk 501 compared with that of the 
elliptical galaxy NGC 5831 observed in both spectral ranges. 

Since early-type galaxies exhibit some gradients in velocity dispersion (Davies et al. 
1983; Fisher, Illingworth & Franx 1995), the measured value of a depends somewhat on the 
distances of the galaxies and the size of the used aperture. In order to compare our values of 
a with the data available in literature (in particular with the MF01 relationship), we applied 
aperture corrections according to the procedure given in Jorgensen, Franx & Kjaergaard 
(1995). The individual measurements of a are therefore corrected to a circular aperture with 
a metric diameter of 1.19/i -1 kpc, equivalent to 3.4" at the distance of the Coma cluster to 
derive central velocity dispersion cr c , given in column 3 of Table 2. When measurements of 
a in two spectral ranges are available, the average is reported. 

No previous systematic study of the stellar velocity dispersion in BL Lacs is available 
although recently Barth, Ho & Sargent (2002) report optical spectroscopy for Mrk 501. 
These authors measured a value of o = 372±18 km s" 1 which differs significantly from ours 
(a = 291±13 km s _1 ). Their measurement of <r, derived from the region 5200 - 5600 A, 
could be reconciled with our value, given their large scatter (81 km s _1 ) fitting the data in 
this range. On the other hand the higher a value obtained from the Ca II triplet lines (8498, 
8542 and 8662 A), that are partly blended with telluric absorptions, appears inconsistent 
with our value within the estimated errors. Applying aperture correction to the value of 
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Barth et al the difference of sigma becomes even larger (by ~ 15 km s 1 ). 

3. Results and discussion 

We have adopted the relationship between Mbh and a c found for nearby early-type 
galaxies that is based on optical spectroscopy (MF01): 

M BH = 1.48±0.24 x 10 8 (a/200) 4 - 65±0 - 48 [M ] (1) 

We assume that this relationship is also valid for AGN (see e.g. MF01) and in particular 
for BL Lacs. This is consistent with our imaging studies of BL Lacs (Falomo & Kotilainen 
1999; Urry et al. 2000; Falomo et al 2000), indicating that all our objects are hosted 
by luminous ellipticals. The derived values of Mbh are reported in column 4 of Table 2, 
where the errors are the composition in quadrature of uncertainties in a and in the MF01 
relationship. Using the Gebhardt et al. (2000) relationship instead of the one by MF01, 
tends to yield slightly lower values of Mbh but does not substantially modify our main 
conclusions. The values of M B h in Table 2 span a factor ~20 from 5 x 10 7 M for PKS 
2201+04 to 9 x 10 8 M for Mrk 501. 

As mentioned above, M B h is also correlated, but with a larger scatter, with the lumi- 
nosity of the bulge of the host galaxy. The host galaxy absolute magnitude M R (uncorrected 
for extinction) and the effective radii R e of the 7 BL Lacs are given in columns 5 and 6 of 
Table 2. M B h was thus calculated following the relationship by McLure & Dunlop (2002): 

log M B H = -0.50+0.05 - 2.91±1.23 [M ] (2) 

The corresponding values of Mbh are given in column 7 of Table 2. For most sources the 
difference of Mbh derived with the two methods is within the estimated uncertainty. The 
average values of M B h for our BL Lacs derived, respectively, from a and the host luminosity 
are: <\ogM BH >a = 8.62 ± 0.23 and <\ogM BH >host = 8.66 ± 0.25. 

In two cases (I Zw 187 and PKS 2201+04) a factor of ~3 difference in Mbh is found. We 
note that for PKS 2201+04 we derive a significantly lower velocity dispersion with respect 
to the rest of the observed sources leading to a low M B h- On the other hand, for this target 
a is well determined (good S/N data and the two spectral ranges giving similar result). 

In our sample there are two LBL type and five HBL type BL Lacs (column 2 of Table 
2). With the caveat that the number of studied objects is very small, we find no significant 
difference of M B h between the two types of BL Lacs. 

The measurements of a combined with the effective radii of the host galaxies can be 
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used to estimate the mass of the hosts through the relationship (Bender et al. 1992) : 

M host =5a 2 r e /G (3) 

This dynamical mass (column 8 of Table 2) turns out to be in the range of 1 - 4 x 
10 11 M Q . The ratio between Mbh and Mhost is in the range of 0.5 - 3.6 x 10~ 3 , with 
< M BH I M host > = 1.4 x 1CT 3 . This is in good agreement with values derived for both AGN 
and inactive galaxies (< M BH /M host > = 1.2 x 1(T 3 ; McLure & Dunlop 2001; MF01). 

In the unified model of radio-loud AGN, BL Lacs are believed to be drawn from the 
population of radio galaxies according to orientation effects (e.g. Urry 1999). It is therefore 
interesting to compare orientation-independent properties, such as the velocity dispersion of 
the host galaxy, of BL Lac and radio galaxy populations. We show in Fig. 2 the comparison 
of the host luminosity vs. loga (the Faber- Jackson relationship) for the BL Lacs with 
respect to a large sample of low redshift radio galaxies (Bettoni et al. 2001). Both samples 
follow quite well the original Faber & Jackson (1976) correlation. The similarity of the 
distributions of a for these two samples (Fig. 3) furthermore implies that the distributions 
of Mbh i n radio galaxies and BL Lacs are indistinguishable, consistent with the model that 
both types of AGN belong to the same population but are observed from different orientation 
angles. 

We thank D. Bettoni for helpful discussions and suggestions on the measurements of 
the velocity dispersion. This work has received partial support under contracts COFIN 
2001/028773 and ASI-IR-35. Nordic Optical Telescope is operated on the island of La Palma 
jointly by Denmark, Finland, Iceland, Norway, and Sweden, in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astrofisica de Canarias. We thank the NOT 
staff for kind hospitality and support. The data presented here have been taken using 
ALFOSC, which is owned by the Instituto de Astrofisica de Andalucia (IAA) and operated 
at the Nordic Optical Telescope under agreement between IAA and the NBIfAFG of the 
Astronomical Observatory of Copenhagen. 
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Fig. 1. — Optical spectra of the BL Lac object Mrk 501 (z = 0.034) and of the nearby 
elliptical galaxy NGC 5831 (z = 0.0055). The spectra are normalized to the continuum and 
plotted in the rest frame. 
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Fig. 2. — The host galaxy stellar velocity dispersion a c vs. the i?-band absolute magnitude 
for the BL Lacs (filled circles) and for low redshift radio galaxies (Bettoni et al. 2001; 
open circles). The solid line indicates the original Faber & Jackson (1976) relationship, 
transformed into the i?-band. 
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Fig. 3. — The distribution of the stellar velocity dispersion a c of low redshift radio galaxies 
(Bettoni et al. 2001; upper panel) compared with that of the BL Lacs studied in this paper 
(lower panel). 
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Table 1. Journal of observations and results. 



Object 


z 


Setup 


Exposure 


S/N 


a 








(sec) 




km s _1 


NGC 5831 


0.0055 


A 


600 


45 


167+5 






B 


600 


55 


185+10 


Mrk 421 


0.031 


B 


3600 


40 


220+10 


Mrk 180 


0.045 


A 


3600 


40 


225+10 


Mrk 501 


0.034 


A 


2400 


40 


265+10 






B 


3600 


55 


280+15 


I Zw 187 


0.055 


B 


3600 


30 


253+15 


3C 371 


0.051 


A 


2400 


40 


255+15 






B 


3600 


35 


265+20 


1ES 1959+65 


0.048 


B 


2400 


18 


180+15 


PKS 2201+04 


0.027 


A 


3600 


30 


148+5 






B 


2400 


50 


153+8 



Table 2. Velocity dispersion and BH masses of BL Lacs 



Object 


LBL/HBL 


<7 C 


\og(M BH ) a 


M R 


R e 


hg(M BH ) bulge 


log(M (w) ) 






km s _1 


[M ] 




kpc 


[M ] 


[M ] 


Mrk 421 


HBL 


236+10 


8.50+0.18 


-23.12 


3.4 


8.65 


11.20 


Mrk 180 


HBL 


244+10 


8.57+0.19 


-22.81 


5.0 


8.50 


11.45 


Mrk 501 


HBL 


291+13 


8.93+0.21 


-23.87 


15 


9.00 


11.59 


I Zw 187 


HBL 


253+15 


8.65+0.18 


-22.22 


4.7 


8.20 


11.39 


3C 371 


LBL 


284+18 


8.88+0.20 


-23.67 


2.9 


8.90 


11.32 


1ES 1959+65 


HBL 


195+15 


8.12+0.13 


-22.48 


6.6 


8.30 


11.27 


PKS 2201+04 


LBL 


160+7 


7.72+0.13 


-22.36 


5.1 


8.27 


11.00 



